Stationary states of trapped spin-orbit-coupled Bose— Einstein condensates 
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We numerically investigate low-energy stationary states of pseudospin-1 Bose-Einstein conden- 
sates in the presence of Rashba-Dresselhaus-type spin-orbit coupling. We show that for experi- 
mentally feasible parameters and strong spin-orbit coupling, the ground state is a square vortex 
lattice irrespective of the nature of the spin-dependent interactions. For weak spin-orbit coupling, 
the lowest-energy state may host a single vortex. Furthermore, we analytically derive constraints 
that explain why certain stationary states do not emerge as ground states. Importantly, we show 
that the distinct stationary states can be observed experimentally by standard time-of-flight spin- 
independent absorption imaging. 



Introduction — Ultracold atomic gases have proven to 
be excellent systems to emulate various phenomena aris- 
ing in condensed matter and high-energy physics. How- 
ever, due to the charge neutrality of the constituent 
atoms, no Lorcntz force acts on these systems in the pres- 
ence of an electromagnetic field. This apparently limits 
the use of these systems in simulating phenomena aris- 
ing from the coupling between a matter field and a gauge 
field. Therefore, methods to generate effective gauge po- 
tentials for ultracold atoms have been actively studied in 
recent years One way to create such an artificial 

gauge field, is to couple the atoms with spatially vary- 
ing laser fields @, [H • This scheme is advantageous since 
it can be used to create both Abelian and non-Abelian 
gauge potentials [l], Q, and the resulting field may be 
controlled and shaped by adjusting the laser beams fioj ]. 

These techniques can be used to create an artificial 
gauge potential that mimics the vector potential arising 
from the spin-orbit coupling (SOC), i.e., the interaction 
that couples the spin and momentum degrees of free- 
dom d [1IHT7I ]. A special case of SOC is the Rashba- 
Dresselhaus (RD) coupling which is actively studied due 
to its relevance in spintronics fTsM20j ] . In RD coupling, 
the vector potential is proportional to the spin- 1/2 oper- 
ator of a particle within a plane. In Ref. [12|, it was pro- 
posed that RD-type SOC could be generalized for spins 
larger than 1/2 in cold atomic gases using the so-called 
TV-pod setup, where N laser beams are used to couple 
N internal atomic ground states to a common auxiliary 
state. A tetrapod setup (N = 4) was suggested to gener- 
ate RD-type coupling where the vector potential is pro- 
portional to a spin-1 operator within a plane. The RD 
coupling has been shown to generate non-trivial struc- 
tures in spin- 1/2, spin-1, and spin-2 Bose-Einstein con- 
densates (BECs) diHH]. 

In the homogeneous RD-coupled BEC, the solutions of 
the single-particle Hamiltonian are plane waves. These 
solutions provide insight into the stationary states of a 
trapped system, which can be approximated by super- 
positions of different number of the degenerate single- 
particle solutions. In Ref. 2lJ, it was stated that the 



ground state of the trapped RD coupled condensate has 
two possible phases in the spin-1 case, namely the plane- 
wave (PW) and the standing-wave (SW) phases. In the 
PW phase the ground state of the condensate can be ap- 
proximated by a single plane wave, whereas in the SW 
phase each spinor component consists of two counter- 
propagating coherent plane waves forming a standing 
wave. In Ref. (25j . symmetry properties of the system 
were used to classify the ground-states of a trapped BEC 
with strong RD coupling. It was predicted that in addi- 
tion to the PW and SW states, also exotic lattice states, 
namely, the triangular-lattice state and the square-lattice 
(SL) state emerge as the ground states of the spin-1 con- 
densate. These states are superpositions of three and 
four plane waves, respectively, and they are invariant un- 
der simultaneous discrete spin and space rotations about 
the z-axis. Such vortex lattices cannot be created by ro- 
tating the condensate as the total angular momentum of 
these states vanishes. The SL state has also been pre- 
dicted to occur in rapidly quenched spin-orbit coupled 
Bose gases (29|. Furthermore, in Ref. 25 1, states preserv- 
ing the combined SO (2) spin-space rotational symmetry 
of the Hamiltonian were found to be ground states in 
some regions of the parameter space. These states can be 
approximated by an infinite number of degenerate single- 
particle solutions. 

In this Rapid Communication, we analyze the energet- 
ics of stationary states arising from the RD-type SOC 



in an optically trapped pseudospin-1 BEC. In Refs. [21 | 
and [25| . exotic ground states were predicted to emerge 
but only a few discrete values of parameters were con- 
sidered. In our analysis, we concentrate on the effect of 
the SOC strength on the energies of the stationary states 
for realistic density-density coupling strengths and de- 
termine the ground state of the condensate in different 
regions of the parameter space. The SL state is found 
to emerge as the ground state for strong SOC, irrespec- 
tive of the spin-dependent interactions. Starting from the 
homogeneous approximation, we analytically derive con- 
straints for the possible ground states and neglect the 
states that are found to be energetically unfavorable. 



2 



These unfavorable states include the triangular lattice, 
described above, for the reasons we will explicate below. 
Furthermore, we propose a robust method to observe the 
exotic stationary states by imaging only the particle den- 
sity in time-of-flight experiments. 

Theory — In the tetrapod setup described in Ref. (l2| . 
four laser beams are used to couple four internal atomic 
ground states to a common auxiliary state. This coupling 
gives rise to three degenerate dark states, that is, zero- 
energy eigenstates of the atom-light Hamiltonian that are 
superpositions of the four atomic ground states. The 
condensed atoms reside in the dark states which play the 
role of internal pscudospin degrees of freedom [11], and 
the order parameter of the system takes the form of a 
three-component spinor denoted by * = 

The stationary states of the system are solved from the 
time-independent Gross-Pitaevskii (GP) equation 



H[*]*(r) = pt>(r) 



(1) 



where the effective Hamiltonian of the system is given by 



2m V i 



aAj +9 + V(r) + c * t (r)*(r). 

(2) 

Here m is the atomic mass, V is the optical trapping 
potential, a is the SOC strength, and Co is the density- 
density coupling constant. The vector A represents the 
non-Abelian gauge potential and the matrix is an ef- 
fective scalar potential. In the tetrapod setup = 
a 2 Til (2m) and the vector potential is of the Rashba- 
Dresselhaus form A = F x e x + F y e y [T^]. With this 
vector potential and assuming that the optical trapping 
potential V is cylindrically symmetric, the Hamiltonian 
in Eq. ([2]) is invariant under the combined spin and spa- 
tial rotation about the z-axis R = e %rrF * +l ' yL z , 

The simple Hamiltonian in Eq. ([2]) yields a good de- 
scription of the system in the case that the four in- 
ternal atomic states coupled to the auxiliary state are 
chosen such that the density-density coupling is roughly 
independent of the internal state. For example, if the 
states are chosen from the F — 2 manifold of 87 Rb, 
the atom-atom interaction Hamiltonian is a sum of the 
density-density, spin-spin, and spin-singlet pairing inter- 
action terms [30( . The magnitudes of the latter two terms 
are of order 1% of the first one [3l|. Hence, we obtain 
the Hamiltonian in Eq. ([2|) by keeping only the density- 
density term which retains its form also in the dark state 
basis. Counting on the possibility that an analogous spin- 
orbit interaction is realized in a genuine spin-1 system, 
we also perform calculations in the presence of a spin- 
spin coupling term H ss = c 2 * t (r)^ r *(r) • T [H, [33j . 
where c 2 is the spin-spin coupling constant and F = 
(F x , F y , F z ) is a vector of spin-1 matrices. We con- 
sider both the ferromagnetic (ca < 0) and the antifcrro- 
magnetic (c2 > 0) interactions. 



Significant understanding of our results can be gained 
by considering the homogeneous RD-coupled BEC. The 
minimum-energy eigenstate of the single-particle Hamil- 
tonian, Hi = (-Ti 2 V 2 + 2ihaA ■ V + 2h 2 a 2 ) /2m, is 
a plane wave ^S>k = $fee lfcr [l2j] . where k lies in 

the xy-plane, $ fc = (l, V2e ie », e 2l(>h ) T /2, and 6 k = 
aictaii(ky/k x ). The corresponding energy is given by 
£k = h 2 (k 2 — 2ak + 2a 2 ) /2m. Let us consider states 
which are superpositions of the degenerate single-particle 
ground states 



E* 

i=i 



(3) 



where B is a normalization constant and f(r) models the 
presence of the trapping potential and the effect of repul- 
sive interactions. Thus, |/(»*)| 2 yields a density profile 
characteristic of trapped condensates. It is straightfor- 
ward to show analytically that, irrespective of /(r), the 
SOC energies of the states with arbitrary n are degener- 
ate, Ssoc = J * (n)t {2ihaA ■ V + 2h 2 a 2 ) ^&r/2m = 
h 2 (— ak + a 2 ) /m. However, some of these superposition 
states give rise to density profiles which are energetically 
unfavorable. Hence, we require that the condensate den- 
sity profile should be smooth, i.e., 



V *< n >(r) //(r; 



0. 



(4) 



This requirement assures that the interference of the 
plane waves does not cause spatial variations in the den- 
sity. We note that if the density-density coupling is very 
weak, the energy contribution coming from the density 
modulations might not be large enough to render the 
state energetically unfavorable. However, for typical ex- 
perimental values of Co, Eq. ([?]) should be fulfilled. Equa- 
tion (|4]) is trivially satisfied for n = 1, but for other 
odd values of n, the condition never holds, implying that 
these states are energetically unfavorable due to den- 
sity modulations. For n = 2 the requirement is satis- 
fied if \9\ — 02 1 = 71", where 0j = arctan(fc^//cj). For 
other even superpositions, the requirement is satisfied 
only if the constants a.j in Eq. ([3]) satisfy the condition 
a*a,j + a*, = with i,j = l,2,,..,n, 6%i(j') = +7r, 
and i ^ f , j ^ i' . We conclude that the energetically 
favorable solutions consist of either a single k or one or 
multiple pairs {k, —k}. In particular, the above con- 
dition implies that the counter-propagating plane waves 
must have equal amplitudes. When n — > oo, the ansatz 
in Eq. ([3]) gives a cylindrically symmetric state and the 
order parameter ^ in the cylindrical coordinates (p, 4>, z) 
reads 



* (p, 0, z) 



*i (p, z)e 1 ^-^ 4 
*o (p, z) e 1 ^ 
(p, z) e l < K+1 ) 



(5) 
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FIG. 1. (Color online) Energies of low-energy stationary states for effectively two-dimensional condensates as functions of 
the SOC strength. Panels (a)-(c) represent the pseudospin case C2 = 0, the antiferromagnetic case (C2/C0 = 0.01), and the 
ferromagnetic case (02/ Co = —0.01), respectively. Here we take So = 200. The various stationary states can be approximated 
by superpositions of degenerate plane-wave solutions of the single-particle Hamiltonian. The plane-wave (PW) state can be 
approximated by a single plane wave, whereas the standing-wave (SW) state corresponds to a superposition of two coherent 
counter-propagating plane waves. Furthermore, the square-lattice (SL) state is formed by two pairs of counter-propagating 
plane waves. The dashed lines correspond to cylindrically symmetric states with different values of ft in Eq. {5J. Here, 0.5 a 2 
has been subtracted from the energies E for clarity. 



describing a single spin vortex. Time reversal symmetry 
implies that for each k / there is a degenerate state 
with — k, and thus we consider only k > 0. 

Numerical results — The stationary states are found 
by solving the time-independent GP equation (P). In 
order to obtain the excited states, we solve the time- 
independent GP equation enforcing symmetries charac- 
teristic of each state. For the cylindrically symmet- 
ric stationary states the form of the order parameter is 
fixed by Eq. ((5|). We carried out most of the compu- 
tations for a pancake-shaped condensate and assumed 
Gaussian profile in the axial direction. A truly three- 
dimensional system was simulated for various values of 
a and the coupling constant C2 and the obtained ground 
states agree with the ones obtained in the effectively two- 
dimensional case. We measure length in units of the 
radial harmonic oscillator length a r = \jh/ (mui r ) and 
energy in units huj r . The dimensionless coupling con- 
stants for the effectively two-dimensional case are given 
by c w \/E-kNcl/ a z , and a = a r a. Here, a is the vacuum 
s-wave scattering length and a z = Ti / (rmuj z ) is the axial 
harmonic oscillator length. The exact form of cq is setup 
dependent. 

Figure [1] shows the energies of various stationary states 
as functions of the SOC strength a. For the pseudospin 
case with C2 = the SW and SL states are nearly de- 
generate for the intermediate values of a, whereas for 
the antiferromagnetic case with C2 > 0, the SW state is 
found to be the ground state of the condensate. For the 
ferromagnetic case with C2 < 0, the PW state emerges for 
small and intermediate values of a. We observe that in 
all three cases, there is a region where the ground state of 
the condensate is a cylindrically symmetric vortex state, 
corresponding to n — > 00 in Eq. ([3]). For the ferromag- 
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FIG. 2. (Color online) (a) The magnetization of the square- 
lattice state. The arrows represent the projection of the 
magnetization to the xy-p\a,ne and the ^-component of the 
magnetization is illustrated with the colormap. Black color 
corresponds to small values and white color to large values, 
(b) The density isosurfaces of the spinor components in the 
square-lattice state. The green color corresponds to l^-i] 2 , 
the black color to l^ol 2 and the red one to | \1> 1 1 2 . The field of 
view is 6.6 x 3.3 x 2.4 and the isosurfaces correspond to the 
density 0.02 N/cfi.. The density isosurfaces are capped with 
density colormap where the lightest color corresponds to the 
largest density. The dimensionless value of the SOC strength 
was taken to be a = 5.0. 



netic case, the ground state hosts a Mermin-Ho vortex 
with k = 1 in Eq. ([5]), whereas for the antiferromagnetic 
case the k = state emerges as the ground state for 
small values of a. In the absence of spin-spin coupling, 
both k = and k = 1 vortices appear as ground states. 
Furthermore, for large SOC strengths the lowest-energy 
state is, irrespective of the value of C2, a square lattice 
(SL) state, which can be approximated by putting n = 4, 
Kl 2 = M 2 - |« 3 1 2 = M 2 , and 9j = jtt/2 in Eq. ©. 
We found two of these lattice states: one in which the ^0 
component has a density minimum in the middle of the 
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FIG. 3. (Color online) Spinor components in the square- 
lattice state. The orientation of the axes in all the panels is 
the same. Panels (a)-(c) represent the amplitudes and panels 
(d)-(f) the phases of 9i, $o, ^-1, respectively. The insets 
in panels (a)-(c) represent the absolute value of the Fourier 
transform for the corresponding spinor component. The hor- 
izontal and the vertical axes on the insets correspond to k x 
and k y , respectively. The values of the coupling constants 
were taken to be So = 200, C2 = 0, and a — 5.0. 

trap and a phase singularity on the z-axis, and another, 
shifted so that vf^ has a density maximum at the center 
and no phase singularity on the z-axis. These states are 
nearly degenerate and their energetics are qualitatively 
the same. Hence, we only present results for the former 
one. The magnetization and the three-dimensional par- 
ticle density for the SL state are shown in Fig. [2] We 
observe two kinds of spin vortices: polar core vortices 
and vortices with a ferromagnetic core polarized alter- 
nately in z and — z directions, corresponding to Mcrmin- 
Ho vortices and antivorticcs. This kind of alternating 
vortex structure cannot be created by rotating the con- 
densate as the total angular momentum of the state van- 
ishes. The amplitude, complex phase, and the Fourier 
transform of each spinor component in the SL state are 
shown in Fig. [3l 

The oscillations in the SL state energy curves in 
Figs. [Ha)-(c) are caused by the increase in the vortex 
density as the SOC strength increases. On the other 
hand, the oscillations in the curves for the SW state and 
the radially symmetric states are caused by the increasing 
number of nodes in the density profiles of the spinor com- 
ponents. We also investigated the effect of the density- 
density coupling on the energetics of the SL state. We 
observed that increasing cq increases the critical value 
of a for which the SL state becomes the ground state of 
the system. Hence, for strongly interacting BECs, strong 
SOC is required to render the SL state the ground state. 

Finally, we demonstrate that the exotic stationary 
states of the spin-orbit-couplcd condensates can be ob- 
served by timc-of-ffight experiments. Since the station- 
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FIG. 4. (Color online) Three-dimensional dynamics of the 
square-lattice state after the removal of the SOC and the op- 
tical trapping in the x- and y-directions. The density is in- 
tegrated in the z-direction. Panels (a)-(c) represent the den- 
sity profile at t — 0, t = l/tj r and t — 2/u) r , respectively. The 
maxima of the colormap in panels (a)-(c) are 0.03, 0.0080 and 
0.0035 in units of N/cft.. The values of the coupling constants 
were taken to be So = 1000, C2 = 0, a — 5.0. 



ary states can be approximated by a single plane wave 
or superpositions of standing waves, they are character- 
ized by either a single k or pairs {fc, —k}. Hence, when 
such a state is left to evolve in time after instantaneously 
removing the SOC and the optical trapping in the x- 
and y-directions, one would expect that the condensate 
density begins to move in the directions specified by the 
fc- vectors. Figure [4] shows the three-dimensional dynam- 
ics for the SL state after the xy-trap and the SOC were 
switched off. We observe that the condensate density 
separates into four equal portions moving in the four di- 
rections dictated by the fc-vectors. We also computed 
the temporal evolution for the SW state and observed 
that it separates into two segments moving in opposite 
directions. For the PW state, the condensate density 
moves in the direction specified by the single k, whereas 
the cylindrically symmetric states of Eq. ([5]) expand ra- 
dially. Thus the states can be observed from the total 
particle density by a standard absorption imaging tech- 
nique, regardless of whether the states are obtained using 
a genuine spin-1 condensate or a pseudospin-1 conden- 
sate. Furthermore, measurement of the separation of the 
interference fringes visible in Fig.[4jb) yields information 
on | As | . 

Conclusions — We have computed the energies of var- 
ious stationary states of trapped Rashba-Dresselhaus 
coupled pseudospin-1 BECs as functions of the spin- 
orbit coupling strength a. Our results indicate that for 
weak spin-orbit coupling, states with a single vortex can 
emerge as ground states. With intermediate values of a, 
the plane-wave and standing-wave states were found to 
be the ground states for the cases C2 < and c 2 > 0, 
respectively. For strong spin-orbit coupling, the exotic 
square lattice is the ground state irrespective of the spin- 
spin coupling strength. This indicates that the emer- 
gence of the SL state as the ground state is not affected 
by the setup-dependent spin-spin coupling term. Impor- 
tantly, we suggested a robust method to observe these 
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states by imaging the total particle density of the con- 
densate in typical time-of-flight experiments. 
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